UHE NEUTRINOS: FUSING GLUONS WITHIN DIFFRACTION CONE 
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Currently available estimates of the gluon-fusion effect in ultra-high energy 
neutrino-nucleon interactions as well as in DIS on protons suffer from uncertainty 
in defining the scattering profile function F(b). Indeed, the area, S, in the impact 
parameter space populated with interacting gluons varies by a factor of 4 — 5 from 
one analysis to another. To get rid of uncertainties we specify the dipole-nucleon 
partial-wave amplitude T(6) which meets the restrictions imposed by both the total 
dipole-nucleon cross section and the small angle elastic scattering amplitude. The 
area S becomes a well defined quantity proportional to the diffraction cone slope. 
We solve numerically the non-linear color dipole BFKL equation and evaluate the 
UHE neutrino-nucleon total cross section. Our finding is that the saturation is a 
rather weak effect, < 25%, up to E v 



1 Introduction 



Practical needs of the neutrino astrophysics [I] inspired many papers on the cross section 
a vN (E v ) for the scattering of Ultra-High Energy (UHE) neutrinos on nucleons and nuclei. 
The question of interest is the interplay of unitarity constraints on a vN (E v ) and the 
evolution of QCD parton densities. The UHE neutrinos probe the gluon density in the 
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target nucleon at very small values of Bjorken x. The BFKL [5] distribution of gluons 
grows fast to smaller x, G(x,Q 2 ) oc x~ A , where, phenomenologically, A w 0.3. Hence, 
the neutrino-nucleon cross sections <7„jv ~ E A violating the Froissart bound. 

The original idea of Ref. [3] developed further in Ref. [4j was that an overlap in trans- 
verse space and recombination of partons leads to a slow down of the growth of the parton 
density and finally to the saturation of parton densities. Quantitative QCD analysis of 
the non-linear effects in terms of the gluon density G(x,Q 2 ) was initiated by Ref. [5]. 
More recently, different derivations of the non-linear BFKL equation were presented in 
[5] . The strength of the saturation effect is usually estimated as 

K~a(x,Q 2 )/S, (1) 

where a stands for some gluon interaction cross section. The radius of the area in the 
impact parameter plane, S — irR 2 , within which the interacting gluons are expected to 
be distributed, varies considerably, from R 2 = 16 GeV -2 in [7j down to R 2 = 3.1 GeV -2 
in [5]. Besides, the area S is assumed to be independent of x. However, because of the 
BFKL diffusion property, under certain conditions, S acquires the Regge contribution 
~aJplog(l/x) 0. 

In this communication we specify the dipole-nucleon partial-wave amplitude which 
meets the restrictions imposed by both the total dipole-nucleon cross section and the 
small angle elastic scattering amplitude. The area S becomes a well defined quantity 
proportional to the diffraction cone slope. We solve numerically the non-linear color 
dipole BFKL equation and evaluate the UHE neutrino-nucleon total cross section. Our 
finding is that the saturation is a rather weak effect, < 25%, up to E v ~ 10 12 GeV. This 
result differs from predictions found in extensive literature on the subject [71 IIP]. 

2 The partial-wave amplitude and diffraction cone slope 

Generalization of the color dipole BFKL approach developed in [11] to the equation for 
diffraction slope B(x,r) proceeds as follows [HJ[T3]. In the impact-parameter represen- 
tation the imaginary part of the elastic dipole-nucleon amplitude reads 



A(Z,r,q) =2 / d 2 bcxpHqk)L(£,r,b) 



(2) 
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and the dipole cross section is cr(£,r) —2j ePbr(£, r, b). The diffraction slope B = 
— 2(d\ogA/dq 2 )\ q 2= for the forward cone is 

r) = \(h 2 ) = - J-. | d 2 b b 2 r(£, r, b), (3) 

where r(£, r, b) is the profile function and b is the impact parameter defined with respect 
to the center of the q-q dipole. In the qqg state, the qg and qg dipoles have the impact 
parameter b + p 2 i/2. Then [12] . 



—g^ = WZri Jdp 1 H(p 1 )-^p 2 )\ 
T(t Pl ,b+ ip 2 ) + T(t p 2 , b + i^) - r(f, r, b) 



(4) 



where £ = log(xo/x) and V'(p) is the radial light cone wave function of the qg dipole with 
the Yukawa screening of infrared gluons 

y/C F a s (Ri) p 
4>{P) = -ifip/iic. (5 

Here 2 are the q-g and q-g separations in the two-dimensional impact parameter plane 
for dipoles generated by the q-q color dipole source, r is the q-q separation and K v [x) is 
the modified Bessel function. At r,p\ >2 -C R c and in the as = const approximation, the 
scaling BFKL equation j2] is obtained. 

The BFKL dipole cross section cr(£, r) sums the Leading-Log(l/a;) multi-gluon pro- 
duction cross sections. Consequently, a realistic boundary condition for the BFKL dy- 
namics is the two-gluon exchange Born amplitude at x — xq = 0.03. The running QCD 
coupling as(Ri) must be taken at the shortest relevant distance Ri = min{r, pi] and in 
the numerical analysis an infrared freezing as{q 2 ) < 0.8 has been imposed. In Ref.|14) it 
was found that incorporation of the asymptotic freedom into the BFKL equation splits 
the cut in the complex j-plane into a series of isolated BFKL-Regge poles. Also, in 9 
it was shown that breaking of scale invariance by a running ag(r) supplemented by the 
finite gluon propagation radius R c , changes the nature of the BFKL pomeron from a 
fixed cut to a series of moving poles with the finite Regge slope a'jp of the pomeron 
trajectory j — a(t) — an>(0) + a'jpt. The preferred choice R c = 0.27 fm gives a'jp = 0.072 
GeV -2 , Aip = aip(0) — 1 = 0.4 and leads to a very good description of the data on the 
proton structure function and the diffraction cone slope |15j at small x. 
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In [T3] the diffraction slope B in the elastic scattering of the dipole r on the nucleon 
(with the gluon-probed radius Rn) was presented in a very symmetric form 



B(H,r) = ^r 2 + ~R 2 N + 2a' 1P Z. 



(6) 
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Figure 1: The non-linear correction to the neutrino-nucleon cross section as a function 
of neutrino energy 



3 Non-linear effects 



Following j6], one can add to the rhs of Eq. ((4]) (assuming iV c 3> 1) the term 

-r(£,p 2 ,b + ip0r(£ )Pl) b + ip 2 ). 

The partial-wave amplitude in the impact parameter space reads 
Then, integrating over b in Eqs. (|4l7p yields 



2B(i,r) 



(J) 



(8) 



on 



= J d 2 Pl m Pl )-^( P2 )\ 2 

x Pi) + <r(€, P2) - o-(£,r) 



Pi)a{t,, P2 ) 
47r(Bi+B 2 ) 



exp 



8(Bi+£ 2 ) 



(9) 
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Shown in Fig. 1 is the ratio 5a /a = 1 — a™ l N /a l uN which quantifies the strength of non- 
linear effects in the vN cross section. Here a™ l N represents the observable uN charged 
current cross section obtained from the numerical solution of Eq.([9|). The cross section 
a l uN derives from the solution of Eq. ([9]) with the non-linear term turned off . The vN 



cross section calculated for 1/R C — 0.75 GeV and R N — 12 GeV is shown in Fig. 2, 
where <j v n is identical to a™ l N . We conclude, that the non-linear effect is rather weak, 
5a /a w 0.25 at E v = 10 12 GeV. This result differs considerably from earlier predictions. 
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Figure 2: The charged current neutrino- nucleon cross section as a function of neutrino 
energy 
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